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Abstract. By molecular dynamics simulation, we study the charge inversion 
phenomenon of a rod macroion in the presence of polyelectrolyte counterions. We 
simulate electrophoresis of the macroion under an applied electric field. When 
both counterions and colons are polyelectrolytes, charge inversion occurs if the 
line charge density of the counterions is larger than that of the colons. For the 
macroion of surface charge density equal to that of the DNA, the reversed mobility 
is realized either with adsorption of the multivalent counterion polyelectrolyte or 
the combination of electrostatics and other mechanisms including the short-range 
attraction potential or the mechanical twining of polyelectrolyte around the rod 
axis. 



1. Introduction 

The charge inversion phenomenon takes place due to strong correlations of a macroion 
with small salt ions in solution, which has recently been studied for the physiochemical 
and bio-engineering systems [IlElEllilS.,6,iJini[Tni[ni[T21[ni[n]. Of particular 
interest is its possibility of facilitating the delivery of genes through negative cell walls 
[ 151116) . In our previous studies by molecular dynamics simulations pillU|. we adopted 
both static and dynamic models for the macroions. Specifically, in the dynamical study 
of electrophoresis with explicit (particle) solvent |l()j . we directly proved the charge 
inversion by measuring the drift of the macroion along the external electric field. The 
net charge of the macroion complex was estimated with the use of the force balance 
Q* ~ v^i, where /i is the electrophoretic mobility and v is the solvent friction. 

The purpose of this paper is to demonstrate the possibility of mobility reversal 
under electrophoresis, as the result of charge inversion, for elongated macroions 
including the DNA. Simulation method and parameters are summarized below. We 
take the system of one macroion, many counterions, colons and neutral particles as 
solvent. The units of length, charge and mass are, a, e, and to, respectively (a ~ 1.4A 
in water and to ~ 40 a.m.u.) The rod is assumed to lie perpendicularly to the applied 
electric field which extends fully across the periodic system. It was shown that a finite- 
length rod rotates and aligns parallel to the electric field due to the rod polarization 
by specific condensation of positive ions at one end toward the electric field and the 
negative ions at the other end I17|. Charge inversion occurs similarly in that case, 
but the discreteness of the surface charge is essential since the counterions need to be 
pinned down on the smooth surface of the rod macroion |17| . 
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The macroion has a radius Rq — 5a, negative charge Qq with — 20e or — 80e, 
and mass 2000to, which is surrounded by the N'^ number of counterions of a positive 
charge Z~^e and the N~ coions of a negative charge ^Z^e. The system is maintained 
in overall charge neutrality, Qq + Z~^e — N~ Z~e — 0. The radii of counterions and 
coions are a"*" and a~, respectively, with the counterion radius being fixed a+ = a, 
and that of neutral particles is a/2. The mass of the coions and counterions is m, and 
that of N^, neutral particles is m/2. Approximately one neutral particle is distributed 
in every volume element (2.1a)'' f» (SA)'' inside the periodic simulation domain of the 
side L = 32a w 45A, excluding the locations already occupied by ions. 

We solve the Newton equations of motion for each particle with the electrostatic 
(Coulombic) and Lennard-Jones potential forces under the uniform applied electric 
field E {E > 0). A large number of neutral particles are used to model the 
viscous solvent of given temperature and to treat the interactions among the finite- 
size macroion, counterions, coions and solvent. The Coulombic forces under the 
periodic boundary conditions |T5j are calculated efficiently with the use of the so- 
called particle-particle-particle-mesh algorithm |19l I20| . with (32)'' spatial meshes, 
the real-space cutoff 10a and the Ewald parameter a w 0.262 (the electric field is 
accurately calculated with the Ewald method) . The volume exclusion effects between 
particles (both charged and non-charged) are treated with the repulsive Lennard-Jones 
potential (j^Lj = 4:e[{A/r,jY^ - (A/r.jf] for r,j = |r, - rj\ < 2^/^, and 0lj = -e 
otherwise. Here is the position vector of the i-th particle, and A is the sum of the 
radii of two interacting particles. We relate e with the temperature by £ = ksT, and 
choose ksT = e^/5ea (we assume spatially homogeneous dielectric constant e). The 
Bjerrum length is thus Xb — e^/eksT — 5a, which is 7 A in water. The equations of 
motion are integrated with the use of the leapfrog method, which is equivalent to the 
Verlet algorithm j21j. The unit of time is r = ay/m/e (w Ips), and we choose the 
integration time step At — 0.0 It. 

The Joule heat produced by the external electric field on ions and transferred 
to neutral particles by collisions is drained by a heat bath for neutral particles at 
the boundaries. The thermal bath which screens hydrodynamic interactions is safely 
adopted in the present study since the hydrodynamic interactions are screened at short 
distances in the electrolyte solvent |22[l23j . We confirmed this fact by good agreement 
of our results with and without the heat bath ^U]. We also note that the effects of 
finite length and rotation of the rod macroion were discussed in [171 . 

2. A Rod Macroion with Polyelectrolytes 

2.1. Polyelectrolyte Counterions and Coions 

We showed previously that the mobility of a rod (cylindrical) macroion reversed its 
sign at zero salt, which was enhanced by addition of small amount of monovalent salt 
if the macroion was strongly charged, arod = IQrodl/^nRrodL > 0.04e/a^ (0.33 C/m^) 
|24|. For reference, the surface charge density of the DNA is 0.19 C/m^. The reversed 
mobility of the rod macroion was more persistent to larger amount of monovalent salt 
than the spherical macroion of the same radius and surface charge density. However, 
the rod macroion with surface charge density equal to that of the DNA was not subject 
to mobility reversal. As will be shown in this paper, polyelectrolyte counterions can 
promote overcharging of the macroion |25| . 

Figure ^ shows the charge inversion of a rod macroion when both the counterions 
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Figure 1. The rod-shaped macroion under polyelectrolyte counterions and 
colons. All chains consist of four monomers; each of 30 counterion polyelectrolyte 
chains carries four unit charges e (red monomers). The coion polyelectrolyte is 
either (left) 10 chains consisting of four unit charges (— e) (green monomers), or 
(right) 20 chains of two unit charges and two neutral monomers (white monomers). 



Table 1. The mobility of the infinite rod macroion Qrod = — 80e under 
polyelectrolyte counterions and colons. Among 30 counterion polyelectrolyte 
chains consisting of four unit charges e, 20 chains are required to neutralize the 
macroion. The line charge density of the coion polyelectrolyte is varied, where 
each chain consists of one, two or four unit charges (— e) and interleaving neutral 
monomers for the rest. Also, the case with monovalent salt, ,. = 80 and 

' ' salt 

^aalt ~ shown in the bottom row. The mobility values are in unit of 

10-2/iOr. 



charges 


X chains 


-4e X 10 


-2e X 20 


-e X 40 


no 


salt 


1.3 


3.2 


3.4 


w/ 


salt 


-4.9 


-0.085 


-0.36 



and colons are polyelectrolytes. For simplicity, all chains have the same length of 
four monomers. All the monomers are charged unity e for 30 chains of counterion 
polyelectrolyte (red monomers). For the colons, there are two cases: (a) 10 chains of 
four unit charges (— e) (green monomers), and (b) 20 chains, each consisting of two unit 
charges (— e) and two neutral monomers (white). When the line charge densities of the 
counterion and coion polyelectrolytes are equal as in Fig^a), the mobility reversal is 
small as shown in Table I (the bin of — 4e x 10 and "no salt" ). This is because the well 
adsorbed coion polyelectrolyte cancel charges of the counterion polyelectrolyte. When 
the line charge density of the coion polyelectrolyte is reduced to a half, adsorption 
of coion polyelectrolyte becomes less as shown in Figdb), and significant revered 
mobility is obtained (the bin of — 2e x 20 and "no salt"). Further reduction of the 
line charge density of coion polyelectrolyte yields nearly the same reversed mobility. 
When the monovalent salt whose charge content is twice as large as that of the coion 
polyelectrolyte is added, the mobility of the macroion is non-reversed (the bottom row 
of Table I). Charge neutralization by the coion polyelectrolyte is good when its line 
charge density is comparable to that of the counterion polyelectrolyte. 
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Figure 2. The mobility of a cylindrical macroion against the ionic strength 
of polyelectrolyte counterions n^j = Z^N~^/L^, where /^or ~ 2.8 X 10~* 
(cm/s)/(V/cm) and O.OOSe/a'^ 0.33 M/1. Unit charges placed along double 
helices provide the macroion surface charge with o-^od ~ 0.02e/a^ (0.17 C/m^ 
~ ijuj^a)- (a-) Each counterion chain consists of three trivalent Z-ions (filled 
circles), (b) the same as (a) except that all counterions are isolated spheres (open 
circles). The colons are larger than the counterions a~ /a~^ = 1.5. The arrow 
shows the isoelectric point. 



2.2. A Rod Macroion with Surface Charge Density Equal to the DNA 

Previous studies of charge-inverted rod macroions (polyelectrolyte) dealt with their 
condensation to a charged surface [7], or the adsorption of metal ions on 
themselves except for a very recent study |^. Here, we treat the charge inversion of 
a cylindrical macroion in the presence of polyelectrolyte counterions, where the surface 
charge density of the macroion is approximately that of the DNA, arod ~ 0.02e/a^. 
The macroion is assumed to be an infinite rod with Qrod ~ — 20e, radius Rrod — 5a 
and mass 2000to, lying perpendicularly to the applied electric field. The surface 
charges are provided by two sets of ten discrete charges (— e) aligned helically at the 
depth a below the surface of the rod; the spacing of unit charges is 4.1a along the 
helix contour. The polymer counterions (polyelectrolyte) consist of the monomers 
either of (i) trivalent ions or (ii) monovalent ions. The colons are monovalent 
with larger radius than the counterions a~ /a^ = 1.5, which is favorable for charge 
inversion 24 . The Bjerrum length is e^/eksT = 5a. The mobility normalization 
is iJor = va/{2\Qrod\/RrodL) = vo/4:TTarod ~ 2.8 X 10~'*(cm/s)/(V/cm), where vq is 
the thermal speed of neutral particles. The Z-ion concentration at rizi ~ 0.005/a'^ is 
approximately 0.33 M/1. 

For the counterion polyelectrolyte with 3e — 3e — 3e monomers, the macroion 
mobility is reversed as denoted by filled circles in Fig|21 The reversed mobility 

1/2 

increases with the ionic strength of Z-ions as njj above the isoelectric point (arrow), as 
predicted by theory [7] . Then, the mobility decreases and gets non-reversed since more 
colons condense on the counterion polyelectrolyte as the number of screening colons 
increases with the ionic strength. This is verified in the radial distribution functions in 
FigO where the profile of the integrated charge for the large ionic strength in Figl^Jb) 
has an overshoot and is less prominent. This overshoot is characteristic of non-reversed 
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Figure 3. The distribution functions of the counterions (open bars) and coions 
(shaded bars) for the rod macroion with double hehx charges whose surface charge 
is equal to that of the DNA. These correspond to the filled circles in Figl2l with 
the ionic strength (a) n^j ~ 0.005a~^ and (b) n^j ~ O.OlSa^'^. 




Figure 4. A rod macroion whose surface charges are provided by discrete unit 
charges along double helices with the surface charge density nearly equal to that of 
the DNA, is shown for (a) polyelectrolyte cations of three trivalent ions Af+^ = 17 
(red spheres), and (b) polyelectrolyte of three monovalent ions = 35 (yellow 
spheres). The former corresponds to the filled circle in Fig|^at n^j ~ 0.005/a^. 
The anions (green spheres) are larger than the cations, a~ /a~^ = 1.5. The external 
electric field points horizontally rightward. 



mobility in spite of the peaked charge profile jH]. At the point where the reversed 
mobility terminates, the charge content of the salt in the surface layer is still less 
than that of the macroion surface charge, ZeCsXo ~ jCfrod- Good adsorption of the 
polyelectrolyte Z-ions is seen in the bird's-eye view plot of Fig^a), where the average 
bond length between the chained Z-ions is 3.2a. This behavior agrees with the result of 
the theory |^. Also, it may be related to the electrophoresis experiment that observed 
reversed mobility for the concentration ratio of [polycation] / [DNA] k, 2, where the 
macroion was a chicken embryo DNA and the polycation was monovalent polymer 
with a long hydrophobic chain {CH2)nCH{C^NHi)+C2H^ |2HI- The differences of the 
experiment from the simulation setting, however, are the monovalence and existence 
of the hydrophobic chain, which will be mentioned in the next subsection. 
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For isolated Z-ions with other conditions fixed as above, the mobility reversal 
does not occur (open circles in FiglSJ. The counterion condensation is weak in 
this case, possibly due to desorption of counterions by thermal agitation. The 
correlation energy between the surface Z-ions, calculated with the Wigner-Seitz cell 
radius Rws = (Zj-Ka^od)^^'^ ~ 6.9a, 



is comparable to thermal energy, which is below the threshold of charge inversion. 
Thus, the polyelectrolyte counterions with a rod macroion is more favorable for the 
mobility reversal than with the isolated counterions for the same ionic strength '17'. 

For a strongly charged rod macroion with Orod ^ 0.08e/a^, the mobility is 
reversed and increases above the isoelectric point. The peak mobility occurs around 
nzi ''^ 0.02a~'^ with /i ^ 0.04^or- The Z-ion correlation energy here is twice large 
compared to the previous case, Z'^e'^ /2eRws ^ Q.^ksT {Rws ~ 3.5a). The reversed 
mobility, i.e. net overcharging, increases in proportion to the macroion bare charge. 
The theory of an infinite-length rod adsorption on a two dimensional plane can 
be applied to the present simulation if one makes the following replacements, the 
DNA — > polyelectrolyte, the planar surface the rod macroion, which predicts net 
overcharging |SIIZ| 



Here, the line charge density of the polyelectrolyte rj ~ e/3.2a predicts 2.6 times more 
overcharging for the strongly charged macroion than the weakly charged macroion. 
Our result agrees well with the theory prediction. 

As the second case, the electrophoresis of the same macroion under all unit-charge 
polyelectrolyte is examined. The short polyelectrolyte whose charge content (chain 
length) is less than lOe is not subject to charge inversion. The electrostatic energy 
between the unit charges with the mutual distance 6 ~ 3a is /eh ~ ^ksT. Although 
the Manning-Onsager condensation of colons [21] on the counterion polyelectrolyte 
is not evident in Fig^fb), the adsorption of such polyelectrolyte is not as strong as 
to cause charge inversion. The radial distribution function and the bird's-eye view 
plot of Fig^b) reveal that about a half of such polyelectrolyte chains are located 
away from the macroion surface, thus the net charge of the macroion complex is 
negative. It is mentioned in passing that the present molecular dynamics simulation 
is not reproducing the charge inversion due to optimal adsorption of unit-charge 
polyelectrolyte through fractionalization |3(Jj. 

2. 3. N on- Electro static Mechanisms for Mobility Reversal 

As mentioned above, the unit-charge polyelectrolyte of reasonable length alone does 
not lead to mobility reversal. Nevertheless, we show below that the mobility reversal 
occurs for such polyelectrolytes if another mechanism is superimposed. For example, 
the polymers with hydrophobic tails feel short-range attraction forces around the DNA. 
Here, this effect is modeled with the attractive Lennard-Jones potential of the depth 
IksT, without truncation of the potential at the distance 2^/^cr. A reversed mobility 
is obtained for the polyelectrolyte counterions of 10 unit charges, /x « O.OlS^or- 
An analytical theory showed reversed mobility for the surfactant having a charged 
head group and hydrophobic tails with strong hydrophobicity x = —Qk^T j25) . Our 
molecular dynamics simulation uses a weak attraction potential but requires rather a 



Z^e^/2€Rws ~ 5.3kBT 



(1) 



Q* {r]/\D)/\v,[ri/arod>^D)- 



(2) 
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Figure 5. A cylindrical macroion with long polyelectrolyte cations of unit charges 
(red) and isolated anions (green) has reversed mobility. The surface charge density 
of the macroion with double helix charges, cr,.a^ ~ 0.02e/a^, is approximately that 
of the DNA. The long chains of polyelectrolyte mechanically wind around the 
rod. The radius and number of anions are a~ = 1.5a and N~ = 60, respectively 
[neutral particles are not shown]. 

long chain. More works are necessary to reconcile these results with the experiment 
using the DNA as the macroion and monovalent counterion polyelectrolyte 

The following simulation result gives an implication where large reversed mobility 
fi ~ 0.073^0 occurs for a very long polyelectrolyte chain of 20e. This mobility reversal 
is due to mechanical twining of the polyelectrolyte chain around the rod axis, as 
depicted in FigEl These chains are not statically attached to the macroion, but they 
pull the macroion by twining toward the electric field on time average. However, when 
the chain length is cut to a half with the total number of counterions fixed, mobility 
is not reversed as the chains can pass by the rod. Although the mechanical twining 
of polyelectrolyte counterions is not due to electrostatics, both the electrostatic and 
non-electrostatic mechanisms can work together to cause the mobility reversal of the 
DNA macroion in the electrophoresis experiments. 

3. Summary 

In this paper, the charge inversion phenomenon of a rod macroion under 
polyelectrolyte counterions was studied. When both the counterions and colons were 
polyelectrolytes, the charge inversion took place if the line charge density of the 
counterions prevailed over that of colons. The mobility was reversed for the rod 
macroion of the surface charge density equal to that of the DNA due to polyelectrolyte 
counterions of multi- valence. 

For short chains of counterion polyelectrolyte consisting of unit charges, the 
electrostatic effect alone was not sufficient to cause mobility reversal. Either the short- 
range attraction due to hydrophobicity or the mechanical twining of the polyelectrolyte 
around the rod macroion could cause the mobility reversal. In real environments, the 
electrostatic mechanism and other effects such as hydrophobic attraction by specific 
configurations or mechanical twining of the polyelectrolytes may collaborate to result 
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in the mobility reversal under the electrophoresis. 
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